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ABSTRACT 

The distribution of warm molecular gas (1000-3000 K), traced by the near-IR H 2 2.12 jim 
line, has been imaged with a resolution < 0.5 arcsec in the central 1 kpc of seven nearby 
Seyfert galaxies. We find that this gas is highly concentrated towards the central 100 pc and 
that its morphology is often symmetrical. Lanes of warm H 2 gas are observed only in three 
cases (NGC 1068, NGC 1386 and Circinus) for which the morphology is much wider and 
extended than the dust filaments. We conclude that there is no one-to-one correlation between 
dust and warm gas. This indicates that, if the dust filaments and lanes of warm gas are radial 
streaming motions of fuelling material, they must represent two different phases of accretion: 
the dust filaments represent a colder phase than the gas close to the nucleus (within ^100 pc). 
We predict that the morphology of the nuclear dust at these scales should resemble that of the 
cold molecular gas (e.g. CO at 10-40 K), as we show for CenA and NGC 1566 by Atacama 
Large Millimeter/submillimeter Array (ALMA) observations, whereas the inner H 2 gas traces 
a much warmer phase of material identified with warmer (40-500 K) molecular gas such as 
CO(6-5) or HCN (as shown by ALMA for NGC 1068 and NGC 1097). We also find that X-ray 
heating is the most likely dominant excitation mechanism of the H 2 gas for most sources. 

Key words: techniques: high angular resolution - astrometry - galaxies: nuclei - galaxies: 
Seyfert - infrared: galaxies. 


1 INTRODUCTION 


Understanding how supermassive black holes (SMBHs) grow and 
how this is connected to the formation and evolution of their host 
galaxy is one of the major questions in astrophysics (e.g. see review 
by Kormendy & Ho 2013). The accretion of gas coming from 
galactic scales (~10 kpc) is considered to be the main fuelling 
mechanism of SMBHs and eventually of active galactic nuclei 
(AGN). In order to trigger and feed nuclear activity, gas must be 
transported from kpc scales to a few Schwarzschild radii around 
the SMBH and lose its large angular momentum. At kpc-scales, 
galaxy interactions and mergers and large-scale bars have been 
proposed as the most efficient mechanisms driving gas into the 
inner kpc region (e.g. Hopkins et al. 2006; Di Matteo et al. 2008; 
see Combes 2006 and Jogee 2006 for a review). On smaller scales 
(a few 100 pc), nested bars (Shlosman et al. 1989), nuclear spirals 


(Maciejewski 2004) and radiation pressure driven inflows (Wada 
2012) are proposed mechanisms for driving matter to the centre. 

A nuclear spiral produced in a gaseous medium by shocks 
induced by the large-scale bar may also drive the gas to the 
central parsecs of the galaxy (e.g. Englmaier & Shlosman 2000; 
Prieto et al. 2005; Storchi-Bergmann et al. 2007; Davies et al. 
2009). Dusty spirals are also a common phenomenon in both 
Seyfert and “normal” galaxies (Martini et al. 2003), and provide 
an alternative mechanism for the fuelling process. The presence 
of nuclear dust filaments connected to kpc-scale dust structures 
was observed in several Seyfert galaxies (e.g. Malkan et al. 1998; 
Hunt & Malkan 2004; Prieto et al. 2014), suggesting that these 
are the channels of inflow of material from the outer parts of the 
galaxy to the nuclear region. Probing a connection between the 
dusty filaments and gas distribution at scales of a few 100 pc 
and AGN feeding has been possible only in a few nearby AGN 
using high-angular resolution observations of the cold molecular 
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gas (e.g. Garcia-Burillo et al. 2009; Krips et al. 2011; Combes 
et al. 2013, 2014) or of warm gas tracers such as the H 2 molecular 
gas (e.g. Muller Sanchez et al. 2006; Neumayer et al. 2007; Davies 
et al. 2009; Hicks et al. 2009, 2013; Mazzalay et al. 2013). 
The latter is usually traced by integral-field spectroscopy (IFU) 
observations limited to a small instrumental field of view (FoV), 
often of less than 3 arcsec x 3 arcsec. For instance, evidence of a 
close relation between dust and molecular gas has been observed 
in the Seyfert galaxies NGC1566 and NGC1433, for which 
recent Atacama Large Millimeter/submillimeter Array (ALMA) 
observations reveal a nuclear spiral of cold CO(3-2) molecular gas 
that is fuelling the nucleus and is spatially coincident with the dust 
structures seen at scales ^ 400 pc (Combes et al. 2013, 2014). In 
NGC 4501, the warm H 2 molecular gas traced by VLT/SINFONI^ 
within a FoV of 3 arcsec x 3 arcsec correlates very well with 
the filamentary dust lanes observed at scales <100 pc (Mazzalay 
et al. 2013). It should be noted though that kinematic studies of 
the molecular warm gas provide controversial results concerning 
the presence of molecular outfiows in the central regions of some 
AGN, as is the case of NGC 1068, for which the results of the 
SINFONI observations were interpreted as an infiow of warm H 2 
gas towards the nucleus at scales of a few pc that could be fuelling 
the AGN (Muller Sanchez et al. 2009) but as an outfiow by Barbosa 
et al. (2014). Unfortunately, the small FoV of IFU observations is 
in many cases insufficient to fully characterize the gas kinematics 
as a large fraction of gas may reside outside the covered IFU FoV, 
as it is shown in this work. 

In this paper we show the full extension of the warm 
(~ 1000-3000 K) H 2 gas of a nearby sample of five type 2 Seyfert 
galaxies (CenA, Circinus, NGC 1068, NGC 1386 and NGC 7582), 
one type 1.5 Seyfert (NGC 1566) and one type 1 LINER/Seyfert 
(NGC 1097), finding that the warm molecular gas is often detached 
from the central dust lanes and filaments that are observed at 
equivalent angular scales as the H 2 gas (~ 0.1-0.4 arcsec). This 
may be pointing to two modes of accretion: a cold mode, traced 
by the dust filaments, and a warm mode, traced by the warm 
H 2 gas and warm molecules such as HCN or CO of high J 
levels traced by ALMA. The existence of two distinct accretion 
modes, hot (10^“® K) and cold, has been shown by observations 
(e.g. Hlavacek-Larrondo et al. 2013), predicted by high-resolution 
simulations of black hole accretion on sub-kpc scales (e.g. Gaspari 
et al. 2013; Bourne et al. 2014) and implemented in cosmological 
simulations by Steinborn et al. (2015) in an attempt to better 
describe a smooth transition between the so-called radio- and 
quasar-modes of accretion. The existence of an intermediate warm 
gas mode of accretion should therefore take place. We tackle 
this aspect in this paper by investigating the association between 
relatively warm (~ 1000-3000 K) H 2 gas and cold dust filaments 
seen in detail in the central few 100 pc of active galaxies. The 
targets have been drawn from a larger sample that includes the 
nearest and brightest AGN that could be observed at high resolution 
in the infrared (IR) using the VLT and adaptive optics (AO; Prieto 
et al. 2010; Reunanen et al. 2010). We present near-IR narrow-band 
H 2 imaging observations with the VLT of this sample of seven 
nearby AGN. The large FoV of the observations (up to 26 arcsec x 
26 arcsec) and high spatial resolution (down to 0.1 arcsec) allows 
us to study the distribution of the molecular H 2 warm gas and its 
potential as a tracer of dust at scales from a few tens of pc to a 


^ Spectrograph for INtegral Field Observations in the Near Infrared 
(SINFONI) on the Very Large Telescope (VLT). 


few 100 pc as well as to estimate the amount of gas available for 
accretion on to the SMBH. In addition to the narrow-band imaging 
data, we present near-IR broad-band imaging data taken with the 
same instrument, optical I- and V-band data from the Hubble Space 
Telescope (HST), AO VLT/SINFONI data. X-ray imaging data 
from the Chandra X-ray satellite, and CO and HCN intensity maps 
from ALMA. 

The paper is organized as follows. The observations and data 
reduction are presented in Section 2. The data analysis is described 
in Section 3. The results obtained are presented and discussed in 
Section 4, while final conclusions are summarized in Section 5. 
An individual description of each target galaxy is provided in 
Appendix A. 


2 OBSERVATIONS AND DATA REDUCTION 
2.1 Adaptive optics IR data 

The seven galaxies listed in Table 1 were observed with the near-IR 
AO assisted instrument Naos-Conica (NaCo; Rousset et al. 2003) 
installed at the Nasmyth focus on ESO VLT (Paranal, Chile). 
Near-IR images in several broad- and narrow-band filters in the 
X-band were taken with NaCo using the pixel scale of 0.027 arcsec 
pixel” A EoV of 26 arcsec x 26 arcsec of the central few 100 
pc was covered by NaCo for all the galaxies. In all cases the AO 
correction was done on the nucleus or a star cluster close to the 
nucleus using wavefront sensing in the optical. Eor the purpose 
of this paper, we make use of the 2 qm broad-band Ks filter 
(Ac = 2.180 qm, AA = 0.350 qm) to obtain the continuum image. 
This should provide the maximum contrast between the nuclear 
emission and the host galaxy (e.g. Prieto et al. 2010, 2014). 

Eor all the targets, narrow-band images were taken with two 
filters, each with different width but both centred on the H 2 1-0 
S(l) 2.12 qm line. The filters used were NB2.12 (width 0.022 qm) 
and IB2.12 (0.06 qm). Exposure time in both filters was the same 
for each galaxy, ranging from 4 to 5 min on-source. Since the 
seven target galaxies are nearby and have small systemic velocities 
(see Table 1), the H 2 1-0 S(l) 2.12 qm line is fully covered by 
the NB2.12 filter. The resolution in these narrow-band images 
was estimated from point-like sources detected in the EoV. In 
NGC 1566 foreground stars were used for resolution estimate, 
while for the rest the estimate was based on stellar clusters, which 
are unresolved with 8-m telescopes in the near-IR. Eor each galaxy, 
the spatial resolution achieved in the two filters is comparable and 
ranges from full width at half maximum (EWHM)~0.12 arcsec for 
NGC 1566 to EWHM-0.35 arcsec for NGC 1097 (see Table 1). 

Eor the galaxies CenA, Circinus and NGC 1386, the 
intermediate-band line-free filter IB2.06 (Ac = 2.060 qm, AA = 
0.060 qm) was also used to estimate the continuum level. 

The NaCo data reduction was performed following 
the standard procedure of sky-subtraction, flat fielding and 
combination of frames with the ESO ECLIPSE package (Devillard 
1999). 

The imaging data are complemented with lEU VLT/SINEONI 
data for NGC 1386, NGC 1566 and NGC 7582. The data reduction 
was performed using the SINEONI pipeline provided by ESO and 
includes: correction for bad pixels, fiat-field, geometric distortions, 
wavelength calibration, reconstruction of the data cube from the 
image spectral slices, background subtraction using sky frames, 
and flux calibration using a standard star. 
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Table 1. Galaxy properties, filters and errors. 


Object 

Galaxy 

AGN 

^bol 

Dl 

Linear 

Ref. 

Dust map 

H 2 gas 

FWHM 

Clusters 

Errors 


type 

type 



scale 


filters 

filters 

Ks 

NB2.12 

registration 

registration 




[ergs“^] 

[Mpe] 

[pc arcsec“^] 







[mas] 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

(12) 

(13) 

CenA 

SO 

Sy2 

6 X 10^2 

3.8 

17 

[1] 

K/F814W 

NB2.12-1B2.12 

0712 

0('13 

4* 

60, 120, 60, 70 

Circinus 

SAb 

Sy2 

8.4 X 10^2 

4.2 

19 

[2] 

K/F814W 

NB2.12-1B2.12 

0716 

0('l6 

3 

40, 30, 30 

NGC 1068 

Sb 

Sy2 

8.7 X 10^^ 

14.1 

70 

[3] 

K/F550M 

NB2.12-1B2.12 

0('l4 

0('l3 

3 

30, 30, 10 

NGC 1097 

SBb 

Ll/Syl 

1.5 X 10^^ 

14.2 

69 

[4] 

K/F814W 

NB2.12-1B2.12 

0714 

0735 

7 

20,30,30,30,20,30,30 

NGC 1386 

SBa 

Sy2 

1.2 X 10^2 

15.3 

74 

[5] 

K/F814W 

NB2.12-1B2.12 

0('09 

0('19 

2* 

10, 10 

NGC 1566 

SABbc 

Syl.5 

5.5 X 10^2 

20.5 

99 

[6] 

K/F438W 

NB2.12-1B2.12 

0('13 

0('l2 

3 

40, 40, 10 

NGC 7582 

SBab 

Sy2 

1.8 X 10^^ 

19.9 

96 

[7] 

F160W/F606W 

SINFONI** 

0716 

0('24 

4 

20, 50, 70, 20 


Column designation: (1) Object name; (2) galaxy type; (3) AGN type: Sy stands for Seyfert, L for LINER; (4) bolometric luminosity, from Mezcua & Prieto 
(2014) and Prieto et al. (2010); (5) luminosity distance; (6) linear scale; (7) references for column (5): [1] Harris et al. (2010); [2] Freeman et al. (1977); [3] 
Marco et al. (1997); [4] Tully et al. (2009); [5] Jensen et al. (2003); [6] NED adopting ifo = 73 km s"! Mpc"!; [7] Terry et al. (2002); (8) and (9) filters 
used to construct the dust maps and H 2 gas maps, respectively; (10) and (11) FWHM of the Ks-hmd and NB2.12 filters; (12) number of stars/clusters used in 
the alignment; (13) error on the position of each star/cluster used in the alignment. 

* for CenA and NGC 1386, the NB2.12 and IB2.12 images have been aligned to the NaCo/^fiT-band image using the nucleus; ** for NGC 7582 no good 
quality map with NaCo is available. The H 2 molecular gas map of this galaxy is thus derived from SINFONI data. 


2.2 HST optical data 

To construct the colour dust maps, we searched in the HST archive 
for the reddest image available for each galaxy. This should provide 
the largest number of common sources to the 2 p.m NaCo images 
to use as reference sources in the image registration. Broad /-band 
images corresponding to the HST/F814W filter from either the 
WFPC2^ or the WFC3^ cameras were retrieved for most of the 
target galaxies (see Table 1). When not available, the V-band 
image from the ACS/F550M filter was used. In the case of 
NGC 7582, the colour maps were built using the WFFC2/F606W 
and NICMOS^/F’/dOW images to avoid the strong contribution of 
the nucleus in the K-hand. In all cases the pixel scale ranges from 
0.04 to 0.05 arcsec pixeP^. 

For data sets with multiple exposures, the retrieved images 
were combined and corrected for geometric distortions using the 
MultiDrizzle package in PyRAF (Koekemoer et al. 2003). Cosmic 
rays removal was applied for those datasets with single exposures 
using the L.A. Cosmic algorithm (van Dokkum 2001) in its 
Python module version^. Further analysis was performed with 
iRAF/PyRAF. 


3 ANALYSIS 

3.1 Extraction of H 2 -line maps 

The extraction of the continuum-free H 2 images was performed 
in two steps. First, the subtraction of the two images centred on 
the 2.12 p.m H 2 line (IB2.12-NB2.12) produced a pure continuum 
(offline) image. This pure-continuum image was further used to 
extract pure H 2 images from the narrow-band NB2.12 image using 
a scaling factor. This process was repeated iteratively, changing the 
scaling factor, until the residuals of the H 2 emission-line image 
were smooth. Effectively this corresponds to i /2 = NB — c x 
{IB — NB), where c is the scaling factor. This process was 


^ Wide-Field Planetary Camera 2. 

^ Wide-Field Camera 3. 

^ Near Infrared Camera and Multi-Object Spectrometer. 

^ http://python.org 

http://obswww.unige.oh/~tewes/cosmics_dot_py/ 


performed to obtain the H 2 maps of all the galaxies in the sample 
except for NGC 7582 (see below). 

For CenA, Circinus and NGC 1386, H 2 maps were also 
extracted using the filters IB2.12 and IB2.06, containing the 
emission line and adjacent continuum, respectively. A similar 
iterative process was applied to extract the H 2 emission-line maps. 
The latter were compared to those based on the NB2.12 filter. In all 
cases, similar H 2 structures are derived using these two approaches. 

In addition to the NaCo narrow-band imaging, we show for 
all sources (Figs. 3-9) the molecular H 2 gas distribution obtained 
from AO VLT/SINFONIIFU observations. For CenA, NGC 1068, 
NGC 1097, NGC 1386 and NGC 1566, the H2 line maps have been 
extracted from the SINFONI data cube using QFITSVIEW®. The 
3D data cube allows an optimized extraction of the line emission, 
subtracting the adjacent continuum measured at both sides of the 
line. The warm SINFONI H 2 maps for Circinus and NGC 1068 
are directly from Muller Sanchez et al. (2006) and Muller Sanchez 
et al. (2009), respectively. 

For NGC 7582, the IB2.06 filter could not be used as reference 
for the continuum level due to contamination from the He I 
emission line at 2.058 p.m, associated with the circumnuclear 
starburst. Thus, the H 2 line map shown in Fig. 9 was extracted from 
the SINFONI data cube following the same procedure as for CenA, 
NGC 1097, NGC 1386 and NGC 1566. 

3.2 Image registration 

For each target galaxy, the HST images were registered to the NaCo 
/f-band image and the H 2 narrow-band images using at least two 
point-like reference sources (usually unresolved stellar clusters in 
the galaxy and avoiding the nucleus when possible) in the common 
FoV to all images. In general, up to three to four reference sources 
could be identified (see Table 1) due to the limited FoV and the 
strong galaxy emission in the near-IR. The image registration was 
performed with the IRAF task imalign and crosschecked with 
an alignment script written by us in IDL^ following the same 
procedure as in Prieto et al. (2014). The relative offsets between the 
images were determined using the centroid position of the reference 

® http://www.mpe.mpg.de/~ott/QFitsView/ 

^ Interactive Data Language. 
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Figure 1. NGC 1566. HST/F438W image with NaCo/ZsT-band continuum 
contours in blue. The FoV is 4 arcsec x 4 arcsec. The colour scale is in 
arbitrary units. North is up and east is to the left. 


sources identified in the common FoV and never using the nucleus 
to register the near-IR and optical images (see Prieto et al. 2014). 
The near-IR nucleus was only used to align the .^f-band image to 
the H 2 narrow-band image in the case of CenA, NGC 1386 and 
SINFONI for NGC 7582. The average of these offsets was used to 
shift and align the images. For each target, the alignment error of 
each reference source is estimated as the standard deviation of its 
position in the different images (Table 1, column 13). The mean of 
these individual errors is taken as the final error of the alignment 
for each target galaxy, which ranges between 20 and 80 mas. 
Therefore, any uncorrected distortions or wavelength-dependent 
effects affecting the centroid position of the point-like reference 
sources used in the image registration are accounted for in the final 
alignment error. See Prieto et al. (2014) for further details. 

The accurate registration of the .^f^-band and HST optical 
images allows us to confirm the shift between the HST/F814W 
and NaCo/^^f^-band peaks previously reported for the type 2 Seyfert 
Circinus and to constrain it to 160 zb 30 mas (Mezcua et al., in 
preparation). For the same galaxy we also note a shift of 40 zb 
30 mas between the peak of the warm H 2 molecular gas and the 
.^f^-band peak. In addition, we report for the first time an offset 
between the optical (HST/F438W) and NaCo/TT^-band peaks of 
emission for the type 1.5 Seyfert NGC 1566 of 50zb30 mas (Fig. 1). 

The NaCo warm molecular H 2 gas emission map of Circinus, 
NGC 1097, NGC 1386 and NGC 7581 has been also aligned to their 
archival broad-band X-ray Chandra image (bottom right panel in 
Figs. 4, 6, 7, and 9). The image registration has been performed by 
matching the position of the H 2 peak to that of the nuclear X-ray 
peak due to the lack of point-like references sources in the common 
FoV of the H 2 and Chandra images. In the case of NGC 1068, 
the SINFONI warm H 2 gas map has been aligned to the Chandra 
High Resolution Camera (HRC) image of the nuclear region of 
NGC 1068 (Wang et al. 2012) by matching the position of the 
nucleus in the X^-band to the nuclear X-ray peak (Fig. 5, middle 
right panel). This same procedure (i.e. matching the position of 
the nucleus) has been performed to align the SINFONI warm H 2 
gas map of NGC 1068 to its ALMA CO(6-5) map (Garcia-Burillo 


et al. 2014) and the NaCo warm molecular H 2 gas emission image 
of NGC 1097, NGC 1566 and CenA to their ALMA HCN (Izumi 
et al. 2013; Martin et al. 2015), CO(3-2) (Combes et al. 2014) and 
CO(2-l) maps (Espada 2013), respectively. 

3.3 Dust maps 

I - K or V - K colour maps of the central few 100 pc of each 
target galaxy were constructed using the ratio between the aligned 
NaCo X-band and HST optical images (Figs. 3-9). These maps 
show the presence of dust in the central few 100 pc of each galaxy. 
The dust morphology is in general filamentary, in some cases more 
disc-like or with irregular shape. In CenA, Circinus, NGC 1068, 
NGC 1386 and NGC 7582 (i.e. those type 2 Seyferts) a filament 
of dust is observed to cross the nucleus, unequivocally identified 
as an outstanding point-like source in the NaCo 2 pm X^-band 
images (e.g. Prieto et al. 2010, 2014; Fernandez-Ontiveros et al. 
2012). Such a nuclear dust lane is not observed in NGC 1566 
(Syl.5) nor NGC 1097 (Syl). In NGC 1097 the dust is observed 
to spiral around the nucleus and to recover the same morphology, 
although with much larger contrast, as that originally discovered 
in the high-angular resolution J-K maps of Prieto et al. (2005). In 
some sources (e.g. NGC 1386) the central dust filaments extend to 
~ 1 kpc and are thus associated with the large-scale dust structure 
of the galaxies (see also Prieto et al. 2014). 

For CenA, NGC 1068, NGC 1386 and NGC 1566 the warm 
molecular H 2 gas emission map from SINFONI has been aligned 
to the dust map by matching the position of the H 2 peak to that 
of the nucleus (or NaCo/K-band peak; e.g. Muller Sanchez et al. 
2009). All maps are centred at the X-band peak and are shown in 
Figs. 3-9. 

3.4 Molecular H 2 warm gas light profile 

To quantify the extent of the nuclear H 2 warm gas distribution 
traced by NaCo, circular or elliptical isophotes were fitted to the 
central regions of the H 2 emission-line images (i.e. we exclude the 
H 2 clumps associated with star formation observed in the spiral 
arms of NGC 1097 or the extended starburst of NGC 7582). The 
extent of these isophotes, taken as three times the background noise 
level (cr) of each H 2 emission-line image, is reported in Table 2 (see 
also Fig. 2). A Sersic function (Sersic 1968) was fitted to the light 
profiles of the warm H 2 gas shown Fig. 2. 


4 RESULTS AND DISCUSSION 

4.1 Morphology of the molecular H 2 warm gas 

Molecular H 2 warm gas emission is seen in the central regions 
of all the galaxies discussed in this paper (see Figs. 3-9). The 
emission is resolved by NaCo in all cases, and its morphology 
is varied, ranging from symmetric emission to showing extended 
filaments. In some cases, the warm H 2 gas is also observed as 
clumps of emission in the outermost spiral arms (e.g. NGC 1097, 
see Fig. 6). The gas is thus observed when it is close enough to the 
nucleus as to be heated by the AGN, though it is equally detected 
when heated by star-forming regions (also visible in our NaCo 
X-band images). A circular aperture was able to reproduce the 
morphology of the molecular H 2 warm gas for all sources except 
for NGC 1386 and NGC 1068, for which an elliptical aperture was 
found to better reproduce the observed morphology. In the case 
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Table 2. Properties of the molecular H 2 gas. 


Object 

(1) 

Radius 

[arcsec] 

(2) 

[pc] 

(3) 

Ln2 

[ergs-^ 

(4) 

HWHM 

[pc] 

(5) 

-^warm 

[Mq] 

(6) 

-^cold 

[M©] 

(7) 

-Vcold 

[cm“^] 

(8) 

T^X —rays 

[ergs-^ 

(9) 

CenA 

2.0 

34 

2.4 X 103’’ 

3 

50 

7 X 10® 

9 X 10® 

1.1 X 10"^® 

Circinus* 

1.4 

26 

3.9 X lO^'^ 

4 

40 

1 X 10’’ 

3 X lO-^ 

4.3 X 10^1 

NGC 1068t 

1.5 

106 

1.7 X 1039 

10 

1000 

5 X 10® 

6 X 10^ 

6.5 X 10"^! 

NGC 1097* 

2.0 

138 

1.3 X 1039 

29 

600 

4 X 10® 

6 X lO-^ 

5.2 X 10^0 

NGC 1386t 

4.3 

321 

1.8 X 1039 

100 

1000 

6 X 10® 

5 X 10"^ 

1.3 X 10"^° 

NGC 1566 

2.2 

220 

4.5 X 1039 

27 

3000 

1 X 109 

6 X 10® 

3.7 X 10^2 

NGC7582*-tt 

0.3 

29 

1.6 X 103® 

14 

100 

5 X 10^ 

1 X 10® 

4.5 X 10"^® 


Column designation: (1) Object name; (2, 3) aperture radius used to extract the flux of the warm H 2 molecular gas, in arcsec and pc respectively; (4) 
luminosity of the warm H 2 molecular gas; (5) HWHM of the Sersic funcion fit to the warm H 2 gas distribution; (6) mass of the warm H 2 molecular gas; (7) 
mass of the total (cold) H 2 molecular gas; (8) density of the total (cold) H 2 molecular gas assuming a gas-to-mass ratio /g = 1% for NGC 1097 and /g = 10% 
for the rest; (9) X-ray luminosity from the highest energy range available in the literature for CenA, Circinus, NGC 1068, NGC 1097, NGC 1566, and 
NGC 7582 (Prieto et al. 2010 and references therein). For NGC 1386, the X-ray luminosity has been derived from the 2-10 keV X-ray flux from lyomoto 
et al. (1997) assuming a luminosity distance of 15.3 Mpc. 

* for Circinus, NGC 1097 and NGC 7582, in addition to the AGN-associated warm gas, there are clumps of H 2 gas associated with star formation; ^ For 
NGC 1068 and NGC 1386, elliptical apertures of ellipticity 0.6 and 0.7 were used, respectively; The H 2 molecular gas properties of NGC 7582 come from 
SINFONI data. 



Figure 2. Distribution of the warm H 2 gas provided by the fit of 
circular isophotes (elliptical for NGC 1068 and NGC 1386) to the NaCo 
emission (SINFONI for NGC 7582). Fitting a Sersic function indicates an 
exponential profile for all sources, with a mean HWHM of 27 pc. 

of Circinus, additional clumpy H 2 warm gas emission is observed 
south and north of the nucleus at a 3cr level. The northern clumps 
seem to be associated with star formation, while the southern ones 
are most likely heated by the nuclear source (see Appendix A). 
Taking these clumps into account, we estimate for Circinus a 
total extent of the nuclear H 2 warm gas of 2 arcsec (38 pc) and 
an average extent of the nuclear H 2 warm gas emission for the 
whole sample of 127 pc. This total extent was defined as the 
maximum size of the 3cr distribution in our maps. The fit of a 
Sersic function to the light profiles of the warm H 2 gas (Fig. 2) 
indicates exponential distributions (w~l) for all sources and that 
the warm H 2 gas emission is concentrated at the centre with an 
average half width at half-maximum (HWHM) of 27 pc. This 
HWHM is in agreement with the average extent of the H 2 warm gas 
emission in Seyfert galaxies when traced by IFU observations (e.g. 
Hicks et al. 2009 find an average HWHM~30 pc for a sample of 


Seyfert galaxies that includes Circinus, NGC 1068 and NGC 1097 
observed by SINFONI with a FoV of up to 3.2 arcsec x 3.2 arcsec). 
When using a larger FoV (e.g. up to 8 arcsec x 8 arcsec; Mazzalay 
et al. 2013; Hicks et al. 2013), the molecular H 2 warm gas is 
observed to extend up to a few 100 pc (e.g. 144 pc for NGC 4536; 
Mazzalay et al. 2013; up to 250 pc for the five Seyfert galaxies of 
Hicks et al. 2013) similarly to what we find for our sample when 
observed by NaCo. Hicks et al. (2013) also study the extent and 
luminosity of a matched-sample of quiescent galaxies, finding that 
these present a fiatter luminosity profile and lower H 2 luminosities 
within a radius of ~100 pc than their sample of Seyferts. The 
sample of seven Seyfert galaxies we study in this paper has an 
average nuclear H 2 luminosity of 1.4 x 10^^ erg s“^ (3.6 x 10^ L©; 
see Table 2), which is one order of magnitude higher than that of the 
sample of quiescent galaxies of Hicks et al. (2013) and of the only 
non-AGN (NGC 3351) included in the sample of Mazzalay et al. 
(2013), and a more centrally concentrated nuclear H 2 distribution 
than that of the quiescent/non-AGN. These results are in agreement 
with those found by Hicks et al. (2013), suggesting that the nuclear 
H 2 molecular gas in the AGN constitutes a gas reservoir of material 
available for fuelling the central SMBH (see also Section 4.5). 


4.2 NaCo versus SINFONI 

The large FoV of the narrow-band NaCo imaging allows us to 
trace the distribution of the warm molecular H 2 gas beyond the 
instrument-limited small nuclear regions traced by IFU instruments 
as SINFONI. The different extent and morphology of the molecular 
H 2 warm gas traced by NaCo compared to that traced by a 
SINFONI FoV of up to ~ 2 arcsec x 2 arcsec can be directly 
observed in Figs. 3-8 for CenA, Circinus, NGC 1097, NGC 1068, 
NGC 1386 and NGC 1566. In all these cases the H 2 warm gas 
in the NaCo maps extends beyond the SINFONI FoV, up to 
~ 6 - 8 arcsec (a few 100 pc) for most galaxies. For some 
sources like Circinus, the gas presents a complex morphology 
very distinct to the relatively symmetric distribution traced by 
SINFONI, while in others like NGC 1097 the extension of the 
H 2 emission seems to have a different orientation for NaCo and 
SINFONI. These morphological differences are not caused by 
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a higher SINFONI spatial resolution, as this is in many cases 
comparable or even lower than that of NaCo (e.g. for NGC 1097 
the SINFONI spatial resolution is FWHMsinfoni=0.25 arcsec, 
Davies et al. 2009, while FWHMNaCo=0.14 arcsec, see Table 1; for 
CenA FWHMsinfoni=0.12 arcsec, Neumayer et al. 2007, while 
FWHMNaCo=0.12 arcsec; for Circinus FWHMsinfoni=0.22 
arcsec, Muller Sanchez et al. 2006, while FWHMNaCo=0.16 
arcsec). The differences between the NaCo and SINFONI H 2 
intensity maps can be ascribed to the technique used to extract 
the H 2 line map from the SINFONI data cube: the continuum 
subtraction from the SINFONI data cube (e.g. in the case of 
NGC 1068; Muller Sanchez et al. 2009) relies on the equivalent 
width of the CO bands. Thus, the derived H 2 emission and 
morphology are quite uncertain in the nucleus, where the CO bands 
are almost diluted by the contribution of the non-stellar emission. 
This can seriously affect the morphology of the inner regions, as 
clearly evidenced for NGC 1068: a nuclear hole is observed in 
the warm H 2 emission-line map of Muller Sanchez et al. (2009) 
(see SINFONI map in Fig. 5), which is not present in the NaCo 
H 2 emission map. We note that such nuclear hole is not present 
either in the warm H 2 emission map of Barbosa et al. (2014), in 
agreement with the NaCo map presented in this work. Therefore, 
the NaCo image does recover more H 2 emission than the SINFONI 
maps shown in Fig. 5. 

In order to check for consistency, we extract the H 2 1-0 
S(l) 2.12 ]im line flux within an aperture in the NaCo H 2 image 
equivalent to that of the SINFONI FoV (see Table 3). This is 
performed for those galaxies for which both NaCo and SINFONI 
molecular H 2 gas fluxes are available (CenA, Circinus, NGC 1097, 
NGC 1386, and NGC 1566). For NGC 1068 such a comparison is 
not viable due to the nuclear hole in the SINFONI H 2 emission-line 
map. Given the higher sensitivity of NaCo, the H 2 1-0 S(l) line 
fluxes measured with this instrument are expected to be higher than 
those obtained with SINFONI. This is indeed the case for Circinus, 
NGC 1097 and NGC 1386, for which the NaCo H 2 1-0 S(l) line 
flux is one to three orders of magnitudes larger than that from 
SINFONI. However, for CenA and NGC 1566 the NaCo H 2 1-0 
S(l) line flux is one to two orders of magnitude smaller than the 
SINFONI flux. No polluting lines or features are observed in the 
region of the SINFONI spectrum close to the H 2 1-0 S(l) line, 
indicating that the H 2 1-0 S(l) line flux difference must be caused 
either by the flux calibration of the SINFONI data cubes or the 
continuum subtraction in the NaCo narrow-band Alter. Calibration 
problems in the SINFONI data cube for the particular case of CenA 
could explain flux differences between NaCo and SINFONI of 
up to a factor 4, since the SINFONI observations were acquired 
under variable conditions (Burtscher, private communication). For 
instance, Israel et al. (1990) determined a well-calibrated H 2 1-0 
S(l) flux in a 6 arcsec aperture of (4.2 zb 0.5) x 10“^^ erg cm“^ 
s“^, which is of the same order as the NaCo flux. 

4.3 Comparison with dust and cold molecular gas 

To characterize whether the molecular H 2 warm gas follows the 
dust and/or is coincident with regions of star formation, the H 2 
warm gas contours are plotted on top of the dust maps of each 
galaxy and on top of the ^-band or /-band maps for those galaxies 
with clumps of star formation. An individual description of each 
galaxy and of its H 2 and dust morphology is provided in the 
Appendix A. In CenA, NGC 1097, NGC 1566 and NGC 7582 the 
nuclear warm H 2 gas presents a diffuse, symmetrical morphology 
that does not correlate with that of the dust. The observed 


Table 3. Comparison of NaCo and SINFONI fluxes 




H 2 1-0 S(l) flux 

Object 

Aperture 

NaCo 

SINFONI 


["] 

[erg cm“^ s“^] 

[erg cm“^ s“^] 

CenA 

3x3 

1.5 X 10-14 

3.4 X 10-13 

Circinus 

0.8X0.8 

1.2 X 10-14 

9.0 X 10-15 » 

NGC 1097 

4x4 

5.9 X 10-14 

9.1 X 10-15 

NGC 1386 

4x4 

7.3 X 10-14 

2.2 X 10-14 

NGC 1566 

3x3 

7.7 X 10-14 

5.8 X 10-13 


* The SINFONI H 2 1-0 S(l) line flux of Circinus is taken from Muller 
Sanchez et al. (2006). 

differences are not caused by a different FWHM between the dust 
maps and H 2 (NB2.12) images (see Table 1), but reflect a real 
different morphology. In Circinus (southern emission), NGC 1068 
and NGC 1386 the morphology of the H 2 gas seems to reproduce 
that of the dust filaments close to the nucleus. However, the 
correspondence between the warm H 2 gas and the dust is not 
one to one. The lanes of warm H 2 gas observed in NGC 1068 
and NGC 1386 are wider than the narrow dust filaments and more 
extended. Ailing gaps in the dust morphology. This is also observed 
in other Seyfert galaxies like NGC 7743 (Davies et al. 2014) and 
suggests that gas and dust trace two different phases at scales 
close to the nuclear region: once dust and gas approach the central 
few 100 pc, the gas is first warmed up by the nucleus (and/or 
in some cases jet shocks; see Section 4.4) and thus presents a 
wider and more extended morphology, while the dust remains 
a colder phase with a temperature well below 2000 K. This 
suggests that the dust morphology should resemble that of the 
cold molecular gas (e.g. CO at 10-40 K). The latter can be tested 
in those galaxies for which high angular resolution millimetre 
observations have been performed using ALMA (CenA, Espada 
2013; NGC 1068, Garcia-Burillo et al. 2014; NGC 1097, Izumi 
et al. 2013; NGC 1566, Combes et al. 2014). In NGC 1566, a dual 
spiral structure of cold CO(3-2) gas emission is observed to emerge 
up to 300 pc from the nuclear region and to indeed coincide with 
the dust extinction (Fig. 8, bottom-left panel; see also Combes 
et al. 2014). Such a spiral structure is however not observed in our 
molecular H 2 warm gas emission map despite the large FoV of 
NaCo (Fig. 8, bottom right panel). A good correlation between the 
cold CO(3-2) gas emission and dust is also observed for the Seyfert 
2 galaxy NGC 1433 (Combes et al. 2013) and for NGC 1068 
(Garcia-Burillo et al. 2014). The suggested two phases of accretion 
in which dust remains colder than the gas when approaching the 
nuclear (few 100 pc) scales is also clearly illustrated in CenA: a 
spatial anticorrelation is observed between the warm H 2 and cold 
CO(2-l) gas, with the H 2 Ailing the central ~2 arcsec region devoid 
of cold gas (Fig. 3, bottom right panel), while the cold CO(2-l) 
gas seems to extend in the same direction as the dust filament that 
comes from a few tens of pc to the nucleus (Fig. 3, bottom-left 
panel). In the case of CenA the warm H 2 gas distribution presents a 
shell-like structure that extends along the direction of the radio jet 
(see Fig. 3, middle-left panel), where no cold CO(2-l) gas emission 
is observed (Fig. 3, bottom panels). This warm H 2 gas morphology 
and the lack of cold gas emission in that direction are thus most 
likely caused by jet bow shocks (see Section 4.4). 

While the cold molecular gas distribution resembles the 
morphology of the dust but not that of the warm molecular gas, as 
predicted by the two phases of accretion scenario and observed in 
all the target sources for which ALMA observations are available. 
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the warm (1000-3000 K) and cold (10-40 K) molecular gas are 
part of the same dynamical structure and should thus show similar 
morphologies (or not be completely spatially decoupled) at ’mixed’ 
or ’intermediate’ temperatures of e.g. 40-500 K. A good illustration 
of this is observed in NGC 1068 and NGC 1097: the morphology 
of the warm H 2 gas closely resembles that of the slightly warm 
molecular gas (T >40 K) such as CO(6-5) or HCN, among others, 
as shown in the ALMA images of NGC 1068 (Fig. 5, bottom 
right panel) and NGC 1097 (Fig. 6, bottom-left panel; see also 
Izumi et al. 2013, Martin et al. 2015), while, as expected from 
the two phases of accretion scenario, the warm H 2 gas appears 
disconnected from the cold dust (e.g. Fig. 6, middle-left panel). 


Table 4. Predicted H 2 lAlpm fluxes from the X-ray models of Maloney 
et al. (1996) for a gas density of 10^ cm 


H 2 1-0 S(l) flux 


Object 

[10^2 cm-2] 

Predicted 
[ergcm“^ s“^] 

Observed 
[ergcm“^ s“^] 

Predicted if f=10% 
[ergcm“^ s“^] 

CenA 

2.4 

1.2 X 10“^^ 

1.6 X 10“^^ 

1.2 X 10“^^ 

Circinus 

1.3 

3.9 X 10“^^ 

1.9 X 10“^^ 

3.9 X 10“^^ 

NGC 1068 

0.4 

4.5 X 10“^^ 

7.1 X 10“^^ 

4.5 X 10“^^ 

NGC 1097 

0.04 

6.0 X 10“^^ 

5.4 X 10“^^ 

6.0 X 10“^^ 

NGC 1386 

0.5 

1.4 X 10“^® 

6.6 X 10“^^ 

1.4 X 10“^^ 

NGC 1566 

0.6 

9.1 X 10“^^ 

8.9 X 10“^^ 

9.1 X 10“^^ 

NGC 7582 

0.9 

2.1 X 10“^® 

3.4 X 10“^® 

2.1 X 10“^^ 


4.4 Source of the warm H 2 molecular gas emission 

To explain the excitation of the warm H 2 gas, we look for 
evidence of spatially associated star formation, shocks produced 
by jet interaction or thermal heating by X-rays. Clumps of 
warm molecular H 2 gas emission coincident with regions of star 
formation in the HST /-band, NaCo/ZT-band or ESO/VISIR mid-IR 
bands are observed in the northern spiral arm of Circinus (Fig. 4, 
see Appendix A), the starburst ring of NGC 1097 (Fig. 6) and 
NGC 7582 (Fig. 9). In these regions, star formation is the heating 
mechanism that excites the H 2 molecules. Apart from these clumps, 
no warm H 2 gas emission is observed beyond scales of a few 100 
pc despite the large FoV (~30 arcsec) of the H 2 intensity maps. 
We thus conclude that the H 2 maps obtained with NaCo show 
the majority of the warm H 2 molecular gas present in the studied 
galaxies. 

No clumps of star formation are observed in the nuclear region 
of any of the sources (except for the northern spiral arm of Circinus 
mentioned above), indicating that the central warm H 2 gas is heated 
by the AGN. In nuclear regions, the H 2 molecules are found to 
be typically excited by thermal processes (e.g. Rodriguez-Ardila 
et al. 2004; Mazzalay et al. 2013; Combes et al. 2014) due to 
either shocks (Hollenbach & McKee 1989) or X-ray illumination 
(Maloney et al. 1996). Most of the galaxies studied here were 
analysed in Reunanen et al. (2003) using long-slit near-IR data. On 
the basis of the H 2 2-lS(l)/l-0S(l) ratio, thermal excitation was 
favoured as the dominant mechanism for the nuclear H 2 observed 
in their sample of 14 Seyfert galaxies. In the present study, where 
we trace the whole distribution of warm H 2 within the central kpc 
region, we can state that no young star formation is seen associated 
with the central but spatially extended H 2 emission in any of the 
galaxies: no evidence for it is seen at the centre in our high angular 
resolution K-band images. UV photons are thus not the cause of 
the source of H 2 excitation, in agreement with the Reunanen et al. 
(2003) conclusion. 

We find, however, a fair spatial correlation between the H 2 
and the X-ray morphology, which is particularly evident in those 
cases where the X-ray and H 2 emission is seen spatially resolved 
at the centre: for NGC 1068 and NGC 1386, the nuclear X-ray 
emission presents an extension and morphology that resembles 
that of the warm H 2 molecular gas (middle-right panel in Fig. 5, 
for NGC 1068, and bottom-right panel Fig. 7, for NGC 1386) and 
is coincident with the extent and morphology of the narrow-line 
region traced by the [OIII] ionized gas (Capetti et al. 1997; Bianchi 
et al. 2006; Wang et al. 2012), suggesting that both the diffuse 
X-ray and the [O III] emission are produced by the same gas 
photoionized by the nuclear continuum (e.g. Iwasawa et al. 2003; 
Wang et al. 2012). The extension of the nuclear 0.3-10 keV 
band X-ray emission is also comparable to the nuclear warm H 2 


* Predicted H 2 l.llpm fluxes assuming a conservative Ailing factor of 
10 %. 


gas distribution in NGC 1097, NGC 7582 and Circinus (Figs. 4, 
6, 9, bottom right panel) indicating that the nuclear warm H 2 
molecular gas is most likely excited by the X-ray emission of 
the photoionized gas, as it was also suggested for e.g. the Seyfert 
galaxies NGC 4151, NGC 5548 and NGC 3227 (Rodriguez-Ardila 
et al. 2004; Storchi-Bergmann et al. 2009). The nuclear X-ray 
emission traced by the Chandra X-ray satellite for CenA (Evans 
& Koratkar 2004) is observed to cover an area larger than the 
extent of the nuclear warm H 2 molecular gas. Eor NGC 1566, the 
X-ray emission is unresolved (Fabbiano et al. 1992) and does not 
allow a direct morphological comparison with the warm H 2 gas 
distribution. 

Hard X-ray photons can easily penetrate the molecular clouds 
before being absorbed and heat the medium to excite the Ipm 
vibration-rotation lines of H 2 . We test this possibility from the 
energy balance point of view by comparing the observed H 2 fluxes 
with those predicted in the models of Maloney et al. (1996). These 
models investigate the effects of X-ray irradiation on molecular 
gas and predict H 2 fluxes for a wide range of X-ray impinging 
luminosities and densities. We follow a similar approach as in 
Rodriguez-Ardila et al. (2004). Accordingly, we determine the 
effective ionization parameter, defined in Maloney et al. (1996) as 
follows: 


C 100 



( 1 ) 


where Lx is the hard X-ray luminosity in units of 10^^ erg s“^. 
He the total gas density in units of 10^ cm“^, d the distance in pc 
of the emitting gas to the AGN, and A^Hatt the attenuating column 
density to the X-rays in units of 10^^ cm“^, and we use Fig. 6 in 
Maloney et al. (1996) assuming ne=10^ cm“^. The X-ray nuclear 
luminosities are taken for most sources from Prieto et al. (2010), 
who compile fluxes at energies above 20 keV. These reflect the 
truly X-ray emerging flux from the nucleus, particularly for the 
Compton thick sources in the sample. For NGC 1386, we use the 
2-10 keV X-ray flux from lyomoto et al. (1997). The average X-ray 
luminosity of the sample is ~ 10^^ erg s“^ (see Table 2). The 
attenuating A^h correspond to the extinction values derived by 
Reunanen et al. (2002) and Prieto et al. (2005, 2010, 2014) for the 
central few 100 pc in these galaxies, which is within the radius the 
H 2 emission is observed. We used the conversion factor A^h= (1.79 
zb 0.03) X 10^^ cm“^ Av (Predehl & Schmitt 1995). The values 
of A^Hatt obtained are typically in the range 10^^-10^^ cm“^ (see 
Table 4). The radius of the nuclear extended H 2 emission is taken 
from Table 2. 

The predicted H 2 flux apparently agrees well with the one 
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measured only for CenA. However, considering the uncertainties 
in the determination of the ionization parameter, which could be 
up to an order of magnitude higher if allowed for lower A^Hatf 
the predicted H 2 fluxes are within the range of observed values 
also for Circinus, NGC 1097 and NGC 1566 (see Table 4). In the 
remaining three, the prediction is an order of magnitude less, even 
if ^Hatt is allowed to get reduced by an order of magnitude -this 
is because the dependence of the H 2 flux on ^ is weak for the 
range of ^ derived for these sources. Also much lower iVn would 
expose easily the core of the H 2 molecular cloud to the X-rays 
and would become fully ionized. The underpredicted H 2 fluxes are 
thus surprising given the spatial association between the H 2 and 
the X-rays. Yet, we know that the interstellar medium in galaxies 
is very clumpy, with molecular gas clouds having average sizes 
in the 10-100 pc range, and Ailing factors of a few percent are 
estimated for the cold dense phase (McKee & Ostriker 1977). A 
way to reconcile the above results is thus to assume a Ailing factor 
in the molecular gas region so that the observed emission arises 
from the most clumpy regions inside the molecular gas region that 
have higher densities; this would also explain the relatively large 
radii at which H 2 emission is detected. Our knowledge of the Ailing 
factor of molecular clouds in the environment of AGN is poor, 
yet assuming Ailing factors (/) in the range of 1% to 10%, in 
line with the values adopted by McKee & Ostriker (1977) in their 
three-component medium of the interstellar medium, would bring 
the results to agreement (Table 4, last column). An exception is 
NGC 1386, where the H 2 emission is the largest among the objects 
in the sample. This source has a jet (Nagar & Wilson 1999) of 
size^lOO pc and power 1.2 x 10^^ erg s“^ (Mezcua & Prieto 2014) 
pointing to the south direction and coinciding in location with the 
southern extension in H 2 , thus shocks induced by the jet might in 
this case be the local source of H 2 excitation. 

Shock excitation of the H 2 molecule caused by jet interactions 
could also take place in CenA and NGC 1097, who have a radio 
jet of total size ~2 kpc and ~0.1 kpc, respectively, and values of 
jet power (6 x 10^^ erg s“^ for CenA and 1.5 x 10^^ erg s“^ 
for NGC 1097; e.g. Mezcua & Prieto 2014 and references therein) 
several orders of magnitude larger than the H 2 luminosity. For 
CenA, this is supported by the well-deflned shell-like structure of 
the warm H 2 molecular gas emission, which extends perpendicular 
to the direction of the radio jet, and by the flnding of a large 
H 2 /PAH (polycyclic aromatic hydrocarbon) ratio, which rules out 
UV excitation in photodissociation regions as the source of the H 2 
line excitation (Ogle et al. 2010). Therefore, although CenA is the 
only source in the sample that is comfortably explained by thermal 
excitation by X-rays (see Table 4, column 4), the bow shock 
morphology of the H 2 gas emission and apparent fair alignment 
between the H 2 gas and the jet direction qualitatively points to 
shocks induced by the jets as a plausible additional mechanism 
for exciting the H 2 gas. In NGC 1068, a recent episode of star 
formation (~30 Myr ago; Storchi-Bergmann et al. 2012) was found 
to occur in the ring-like structure of warm H 2 gas located at ~ 100 
pc from the nucleus. According to simulations, the formation of 
such rings and star formation at radii of ~ 100 pc could result 
from the shock waves created by a young radio jet (Gaibler et al. 
2012). The excitation of warm H 2 molecular gas ring-like structure 
in NGC 1068 could thus be also explained by a shock jet interaction 
(see Fig. 5, middle-left panel). In this case, the off-centred nuclear 
position of the ring could be explained by inhomogeneities in the 
gas that surrounds the nucleus (Barbosa et al. 2014). 

The large extent of the nuclear warm H 2 gas, up to scales 
of a few 100 pc, can be thus explained for all the targets 


by thermal processes, in agreement with previous studies (e.g. 
Rodriguez-Ardila et al. 2004; Mazzalay et al. 2013; Combes et al. 
2014). Given the similar symmetrical morphology between the 
nuclear X-ray and warm H 2 molecular gas emission. X-ray heating 
is the most likely excitation source of the nuclear warm H 2 gas for 
most sources, possibly also combined with shocks from radio jets 
for CenA, NGC 1097 and NGC 1068. 

Finally, we note that the nuclear X-ray luminosity of the target 
sources is one to four orders of magnitude larger than the nuclear 
H 2 luminosity (see Table 2), indicating that a small fraction of 
X-ray photons penetrating and heating the central H 2 molecular gas 
cloud is enough to excite the H 2 gas. The X-ray luminosities of the 
AGN studied here are also two to four orders of magnitude higher 
than those of the quiescent galaxies of Hicks et al. (2013) (see also 
Mazzalay et al. 2013), which show not so centrally concentrated 
nuclear warm H 2 gas distributions and lower H 2 luminosities 
than the AGN, thus reinforcing the AGN X-ray emission as the 
dominant excitation mechanism of the warm H 2 molecular gas. 


4.5 Molecular gas mass and density 


The large FoV of the NaCo observations reveals the extent and 
morphology of all the warm molecular H 2 gas present in the studied 
galaxies at scales of a few 100 pc. Following Reunanen et al. 
(2002), the mass of this warm H 2 gas (Mwarm) is quantifled using 
the flux extracted from the H 2 intensity maps and assuming a gas 
temperature of 2000 K, the H 2 l-OS(l) transition probability As(i) 
= 3.47 X 10“^ s“^ (Turner et al. 1977) and the population fraction 
of the (z/, J) = (1,3) level /(i, 3 ) = 0.0122 (Scoville et al. 1982) as: 


Mw 


5.0875 X 

^Mpc 


2 -^1-05(1) ^q0.4277A2.2 

erg s^cm“2 


( 2 ) 


where A 2.2 is the extinction at 2.2 /xm (Scoville et al. 1982). The H 2 
flux was extracted in a nuclear circular aperture for CenA, Circinus, 
NGC 1068, NGC 1097, NGC 1566 and NGC 7582, and an elliptical 
aperture for NGC 1068 and NGC 1386 with a size corresponding to 
the extent of the warm H 2 gas emission given by the isophote fitting 
(Table 2). The values of A 2 . 2 , derived from regions surrounding 
the nucleus, are taken from Reunanen et al. (2002) for CenA and 
NGC 1566, Prieto et al. (2005) for NGC 1097, Prieto et al. (2010) 
for Circinus and NGC 1068 and Prieto et al. (2014) for NGC 1386 
and NGC 7582, and range between A 2.2 = 0.2 and 1.6 for all 

sources. The masses of the nuclear warm H 2 gas are reported in 
Table 2 and range between 40 M© (for Circinus) and 3000 M© 
(for NGC 1566). 

In Section 4.1 we found that the nuclear H 2 luminosity of the 
sample is on average 10 times higher than that of the quiescent 
sources from Hicks et al. (2013) and Mazzalay et al. (2013). The 
mass of the AGN molecular gas reservoir is thus ~ 10 times higher 
for the Seyfert galaxies studied here than for the quiescent galaxies 
(as found by Hicks et al. 2013). This suggests that either these AGN 
have more material available for accretion on to the central SMBH 
than the non-AGN galaxies, which could also explain the higher 
X-ray luminosities of the Seyfert galaxies, or that the gas is more 
effectively heated by the AGN in the Seyfert galaxies, so that in 
these we see the H 2 emission while in the normal galaxies the H 2 
gas remains invisible. 

The total content of molecular cold gas can be estimated 
using a correlation between H 2 2.12 /xm luminosity and cold gas 
mass derived from CO observations (e.g. Dale et al. 2005; Muller 
Sanchez et al. 2006; Mazzalay et al. 2013). This allows us to 
obtain an estimate of the amount of cold gas present in the nuclear 
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regions and thus available for accretion on to the AGN. Using 
the Mcoid-^H 2 relation of Mazzalay et al. (2013), equation 4, we 
obtain that the cold H 2 gas mass ranges between 7x10® M© 
(for CenA) and 1 x 10® M© (for NGC 1566; see Table 2). The 
value of Mcoid obtained for NGC 1566 is in agreement with 
the total molecular mass of 1.3 x 10® M© obtained from CO 
Swedish-ESO Submillimetre Telescope (SEST) observations but 
one order of magnitude lower than the Mcoid obtained from 
smaller EoV CO observations (e.g. 18 arcsec with ALMA; see 
Combes et al. 2014 for a discussion), which indicates that NaCo 
recovers the total content of molecular gas for NGC 1566. Eor 
CenA, Circinus and NGC 1068, the estimated values of Mcoid 
are one order of magnitude lower than those obtained from CO 
observations (Curran et al. 1998; Schinnerer et al. 2000; Espada 
et al. 2009), which can be mainly attributed to the larger EoV of the 
CO observations. We note though that better agreement would not 
be achieved even if we considered a larger aperture for the warm 
H 2 gas, as (except for the clumps in the spiral arm of NGC 1097) 
most of the detected warm H 2 gas is concentrated in the nuclear 
regions and no emission is observed further out despite the large 
EoV ofNaCo. 

Given that the morphology of the nuclear warm H 2 molecular 
gas appears diffuse and symmetrical in most sources, we can 
estimate its density assuming a spherical volume with a gas mass 
fraction (or gas-to-mass ratio) fg=l% for NGC 1097 and /g=10% 
for the rest of the sources (following Hicks et al. 2009). This 
/g represents the fraction of the volume occupied by a certain 
gas phase and should not be confused with the filling factor (/) 
introduced in Section 4.4®. The volume is derived using the same 
radius as that considered when calculating the nuclear H 2 gas 
masses. We obtain cold gas densities ranging Nco\d= 6 x 10® — 
— lx 10® cm“®, which is consistent with the iVH 2 = 10^-10® 
cm“® typically derived from CO observations (e.g. Mauersberger 
& Henkel 1989; Bayet et al. 2009). 


5 CONCLUSIONS 

In this paper we have presented the morphology and properties of 
the warm molecular H 2 gas of a sample of seven nearby Seyfert 
galaxies observed with VLT AO and narrow-band imaging with 
NaCo. Dust colour maps of these galaxies are also constructed 
using the ratio between broad-band ^-band images and V- and 
/-band HST images. The large EoV of the observations has allowed 
us to trace the distribution of the warm H 2 gas and dust up to scales 
of a few 100 pc with a resolution <0.5 arcsec. 

Eor all galaxies, the warm H 2 gas is concentrated in the 
nuclear regions (with an HWHM of 27 pc) and does not extend 
beyond ~ 300 pc from the nucleus. The warm H 2 gas is thus 
visible when the gas is close enough to the heating source (the 
nucleus, where the peak of H 2 emission is observed) and the H 2 
molecules are excited enough to produce ro-vibrational emission 
lines. Eurther out, the H 2 gas is not warm enough and only the cold 
gas is visible through e.g. CO observations. In NGC 1097, Circinus 
and NGC 7582 some clumps or lanes of H 2 emission are in addition 
observed to coincide with regions of star formation. 

The excitation mechanism for all the target sources is found 

® Hicks et al. (2009) assumed a uniform gas distribution in their work (this 
is, a filling factor of 100%) though the authors conclude that the gas must 
be clumpy. This implies a filling factor / <1, in line with the hypothesis 
made in Section 4.4 of assuming a filling factor for the H 2 molecular gas. 


to be thermal. The cospatial location of H 2 and X-rays, also 
coinciding on physical extension in those cases where spatial 
resolution allows for size measurements, is a strong indication for 
nuclear X-rays being the major source of excitation of H 2 . We 
conclude that irradiation by X-rays is the predominant source of 
excitation for most sources under reasonable assumptions of the 
H 2 clouds volume filling factor (of 10%). 

In NGC 1068, NGC 1386, NGC 1566 and Circinus, the 
morphology of the warm H 2 gas traced either by NaCo or 
lEU observations resembles within a few tens of pc that of the 
filamentary dust lanes that extend from the central parsecs to ~ 
1 kpc scale. These dust filaments are potential fuelling channels 
that transport the gas to the nuclear region. This suggests that 
the nuclear warm H 2 gas of these galaxies comes from outer 
(kpc-scale) regions following these same channels of infiow and 
warms up as it approaches the nucleus. The observed lanes of warm 
H 2 gas are however much wider and extended than the filaments of 
dust, which suggests that gas and dust trace two different states of 
accretion with different temperatures. In such scenario the nuclear 
dust morphology at scales of a few 100 pc should resemble that 
of the cold (10-40 K) molecular gas, as we observe for CenA 
and NGC 1566 by ALMA CO observations, while the inner warm 
(1000-3000 K) H 2 gas presents a wider morphology that should 
resemble that of the slightly warmed gas (40-500 K), as indicated 
by warmer molecular tracers (e.g. CO(6-5), HCN) observed in 
NGC 1068 and NGC 1097 with ALMA. 
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Figure 3. CenA. Top left: HST/F814W image with NaCo/^-band continuum contours in white. The position of the nucleus and its error is marked with a 
cross. Top right: NaCo/ZsT-band image. The contours of the warm H 2 molecular gas obtained from NaCo are shown in white. They start at three times the 
background noise level and increase in factors of 2. Middle left: K/F814W ratio or dust map. The contours of the warm H 2 molecular gas obtained from NaCo 
are shown in white. The orientation of the radio jet (e.g. Clarke et al. 1992) is indicated with a dashed blue line. Middle right: Intensity map of the warm H 2 
molecular gas obtained from IFU SINFONI. The contours of the warm H 2 molecular gas obtained from NaCo are shown in black. Bottom left: K/F814W 
ratio or dust map. The contours of the cold CO(2-l) molecular gas obtained from ALMA (Espada 2013) are shown in white. Contours start at three times 
the rms noise of 2.5 mJy beam“^ (Espada 2013). Bottom right: ALMA CO(2-l) integrated intensity map (from Espada 2013). The contours of the warm 
H 2 molecular gas obtained from NaCo are shown in black. The colour scale is linear and starts at three times the rms noise of 2.5 mJy beam“^. The colour 
palettes in the other panels are in arbitrary units. The EoV for all panels is 8 arcsec x 6 arcsec. North is up and east is to the left. 
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Figure 4. Circinus. Top left: HST/F814W image with NaCo/.Sr-band continuum contours in blue. The FoV is 16 arcsec x 16 arcsec. Top right: NaCo/ZsT-band 
image with the same FoV as the previous panel. The contours of the warm H 2 molecular gas obtained from NaCo are shown in blue. They start at three times 
the background noise level and increase in factors of 3. Middle left: K/F814W ratio or dust map with a FoV of 8 arcsec x 8 arcsec. The contours of the warm 
H 2 molecular gas obtained from NaCo are shown in blue. The orientation of the narrow line region (e.g. Marconi et al. 1994) is indicated with a dashed green 
line. Middle right: Intensity map of the H 2 molecular gas obtained from IFU SINFONI (from Muller Sanchez et al. 2006). Bottom left: HST/F814W image 
with a FoV of 8 arcsec x 8 arcsec. The contours of the warm H 2 molecular gas obtained from NaCo are shown in blue. Bottom right: Chandra X-ray image 
in the 0.5-7 keV band with the same FoV as the previous panel. The contours of the warm H 2 molecular gas obtained from NaCo are shown in blue, colour 
scales are linear and in arbitrary units. North is up and east is to the left. 
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Figure 5. NGC 1068. Top left: HST/F550M image with NaCo/ZsT-band continuum contours in blue. The FoV is 8 arcsec x 8 arcsec. Top right: Intensity map 
of the warm H 2 molecular gas obtained from IFU SINFONI (from Muller Sanchez et al. 2009). The FoV is the same as the previous panel. The contours of the 
warm H 2 molecular gas obtained from NaCo are shown in blue. They start at three times the background noise level and increase in factors of 3. Middle left: 
K/F550M ratio or dust map with a FoV of 4.4 arcsec x 4.4 arcsec. The contours of the warm H 2 molecular gas obtained from IFU SINFONI are shown in blue. 
The orientation of the radio jet (e.g. Gallimore et al. 1996) is indicated with a dashed green line. Middle right: Chandra HRC image from Wang et al. (2012) 
with the same FoV as the previous panel. The contours of the warm H 2 molecular gas obtained from SINFONI are shown in blue. Bottom left: K/F550M ratio 
or dust map with a FoV of 8 arcsec x 8 arcsec. The contours of the CO(6-5) molecular gas obtained from ALMA (Garcia-Burillo et al. 2014) are shown in 
blue. Contours levels go from 5a to 240a, where cr = 2 Jy beam“^ km s“^ (Garcia-Burillo et al. 2014). Bottom right: ALMA CO(6-5) integrated intensity 
map (from Garcia-Burillo et al. 2014) with a FoV of 4.4 arcsec x 4.4 arcsec. The contours of the warm H 2 molecular gas obtained from SINFONI are shown 
in blue. The colour scale is logarithmic in 0.87 Jy beam“^ km s“^. The colour palettes in the other panels are logarithmic and in arbitrary units. North is up 
and east is to the left. 
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Figure 6. NGC 1097. Top left: HST/F814W image with NaCo/Z^T-band continuum contours in blue. The FoV is 26 arcsec x 26 arcsec. Top right: NaCo//^-band 
image with the same FoV as the previous panel. The contours of the warm H 2 molecular gas obtained from NaCo are shown in blue. They start at three times 
the background noise level and increase in factors of 2. Middle left: K/F814W ratio or dust map with a FoV of 4 arcsec x 4 arcsec. The contours of the warm 
H 2 molecular gas obtained from NaCo are shown in blue. Middle right: Intensity map of the warm H 2 molecular gas obtained from IFU SINFONI. The FoV 
is the same as the previous panel. The contours of the warm H 2 molecular gas obtained from NaCo are shown in black. Bottom left: HCN image with a FoV 
of 26 arcsec x 26 arcsec. The contours of the warm H 2 molecular gas obtained from NaCo are shown in blue. Bottom right: Chandra X-ray image in the 
0.2-10 keV band with the same FoV as the previous panel. The contours of the warm H 2 molecular gas obtained from NaCo are shown in blue, colour scales 
are logarithmic and in arbitrary units. North is up and east is to the left. 
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Figure 7. NGC 1386. Top left: HST/F814W image with NaCo/ZsT-band continuum contours in blue. Top right:NaCo/^-band image. The contours of the warm 
H 2 molecular gas obtained from NaCo are shown in blue. They start at three times the background noise level and increase in factors of 2. Middle left: 
K/F814W ratio or dust map. The contours of the warm H 2 molecular gas obtained from NaCo are shown in blue. The orientation of the radio jet (e.g. Nagar 
& Wilson 1999) is indicated with a dashed green line. Middle right: Intensity map of the warm H 2 molecular gas obtained from IFU SINFONI. The contours 
of the warm H 2 molecular gas obtained from NaCo are shown in blue. Bottom left: HST/F547M image. The contours of the warm H 2 molecular gas obtained 
from NaCo are shown in blue. Bottom right: Chandra X-ray image in the 0.2-10 keV band. The contours of the warm H 2 molecular gas obtained from NaCo 
are shown in blue, colour scales are logarithmic and in arbitrary units. The FoV for all panels is 10 arcsec x 10 arcsec. North is up and east is to the left. 
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Figure 8. NGC1566. Top left: HST/F438W image with NaCo/ZsT-band continuum contours in white. The FoV is 10 arcsec x 10 arcsec. Top 
right:NaCo//sr-band image with the same FoV as the previous panel. The contours of the warm H 2 molecular gas obtained from NaCo are shown in white. 
They start at three times the background noise level and increase in factors of 2. Middle left; K/F438W ratio or dust map with a FoV of 6 arcsec x 6 arcsec. 
The contours of the warm H 2 molecular gas obtained from NaCo are shown in white. Middle right: Intensity map of the warm H 2 molecular gas obtained 
from IFU SINFONI. The FoV is the same as the previous panel. The contours of the warm H 2 molecular gas obtained from NaCo are shown in black. The 
orientation of the narrow line region (e.g. Schmitt & Kinney 1996) is indicated with a dashed green line. Bottom left: K/F438W ratio or dust map with a FoV 
of 10 arcsec x 10 arcsec. The contours of the cold CO(3-2) molecular gas obtained from ALMA (Combes et al. 2014) are shown in white. Contours levels are 
1.65, 2.7, 6.5, 10, 14 and 17 times 0.87 Jy beam“^ km s“^. Bottom right: ALMA CO(3-2) integrated intensity map (from Combes et al. 2014) with the same 
FoV as the previous panel. The contours of the warm H 2 molecular gas obtained from NaCo are shown in blue. The colour scale is linear in 0.87 Jy beam“^ 
km s“^. The colour palettes in the other panels are in arbitrary units. North is up and east is to the left. 
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Figure 9. NGC 7582. Top left: HST/F160W image with NaCo/ZsT-band continuum contours in blue. The FoV is 5 arcsec x 5 arcsec. Top rightiNaCo/ZsT-band 
image with the same FoV as the previous panel. The contours of the warm H 2 molecular gas obtained from IFU SINFONI are shown in blue. They start at 
three times the background noise level and increase in factors of 1.2. Middle left: HST/F160W image with a FoV of 6 arcsec x 6 arcsec. The contours of the 
warm H 2 molecular gas obtained from IFU SINFONI are shown in blue. Middle right; F160W/F606W ratio or dust map with the same FoV as the previous 
panel. The contours of the warm H 2 molecular gas obtained from IFU SINFONI are shown in blue. The orientation of the narrow line region (e.g. Ulvestad & 
Wilson 1984; Morganti et al. 1999) is indicated with a dashed green line. Bottom left: VISIR/PAH2J2 image (11.88 p.m) with the same FoV as the previous 
panels. The contours of the warm H 2 molecular gas obtained from IFU SINFONI are shown in blue. Bottom right: Chandra X-ray image in the 0.3-10 keV 
band with the same FoV as the previous panels. The contours of the warm H 2 molecular gas obtained from IFU SINFONI are shown in blue, colour scales are 
logarithmic and in arbitrary units. North is up and east is to the left. 
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APPENDIX A: DISCUSSION ON INDIVIDUAL OBJECTS 
CenA 

Centaums A (also NGC5128 or CenA), at a distance of 3.8 zb 
0.1 Mpc (Harris et al. 2010), is the most nearby elliptical galaxy 
and a merger archetype. It hosts a Sy2 nucleus, covered by a 
prominent dust lane (see Fig. 3, top left), that emits a powerful 
radio and X-ray jet up to kpc scales. The NaCo narrow-band 
imaging of CenA shows that the molecular H 2 gas is concentrated 
in a spherical nuclear region of radius ~0.3 arcsec (~5 pc; 
Fig. 3, top right panel). Further out, fainter shells of gas emission 
are traced up to ~3 arcsec (~50 pc) north-east and south of 
the nucleus. This extended morphology is in agreement with 
the overall distribution observed in the SINFONI map (Fig. 3, 
middle-right panel), where shell features are observed (Neumayer 
et al. 2007). These were interpreted as jet bow shocks, which is 
consistent with the orientation of the radio jet at pc scales (Fig. 3, 
middle-left panel). The warm H 2 gas distribution does not follow 
the filamentary dust morphology traced by the — / colour map 
nor the cold CO(2-l) gas emission (Espada et al. 2009; Espada 
2013; Fig. 3, middle and bottom left panels). The nucleus and jet 
of CenA are the most likely heating sources that excite the H 2 
molecules (Bicknell et al. 2013). 


Circinus 

At a distance of 4.2 zb 0.8 Mpc (Freeman et al. 1977), Circinus is 
one of the nearest Seyfert galaxies. It has a one-sided ionization 
cone north-west of the nucleus (Marconi et al. 1994; see also 
Fig. 4, middle-left panel). The nucleus of Circinus, which has 
been resolved by near-IR imaging (Prieto et al. 2004) and by 
interferometric observations at mid-IR (Tristram et al. 2014), is 
obscured at optical wavelengths (Prieto et al. 2004; Mezcua et al. 
in preparation). 

The warm H 2 molecular gas emission is concentrated within a 
region of ~2 arcsec (~38 pc) radius, with the nuclear morphology 
in north-east-south-west direction (perpendicular to the ionization 
cone). This nuclear elongation is also observed in the small FoV 
of SINFONI presented by Muller Sanchez et al. (2006); see also 
Fig. 4, middle-right panel. The NaCo H 2 narrow-band image, with 
a larger FoV than SINFONI, reveals that extended H 2 emission 
in the form of ~1 arcsec-wide arches is also detected ~2 arcsec 
north-east and south-south-west of the nucleus. The northern 
spiral-like structure in the north-east curling turning towards west 
seems to be coincident with clumps of star formation in the 
/-band and Ha image (Fig. 4, bottom-left panel; Mezcua et al. in 
preparation), which suggests that the molecular gas in this region 
is heated by young stars, while the southern clumps seem to follow 
the spiral morphology of the dust colour map. 

NGC 1068 

The spiral SB galaxy NGC 1068, at a distance of 14.1 Mpc (Marco 
et al. 1997), hosts a very bright Sy2 nucleus that is obscured 
from optical wavelengths up to l/xm (e.g. Capetti et al. 1997; 
Prieto et al. 2010, 2014). The warm molecular H 2 gas distribution 
traced by NaCo covers a wider FoV than the SINFONI map and 
recovers the nuclear H 2 emission missing in the SINFONI map 
of Muller Sanchez et al. (2009) (nuclear hole in Fig. 5, upper 
right panel). It also shows clear diffraction effects (i.e. detector 
spikes) caused by the bright nucleus of this galaxy. Despite these 


diffraction effects and the nuclear differences, the extended NaCo 
morphology is consistent with that traced by SINFONI (e.g. Muller 
Sanchez et al. 2009, Barbosa et al. 2014; see also Galliano & Alloin 
2002 for a VLT/ISAAC map): two prominent gas concentrations 
east and west of the nucleus and a collar-like structure south-west 
of it are observed on the NaCO map, similar to the SINFONI 
morphology and to the morphology of the slightly warm molecular 
gas (T ~40-50 K) traced by recent CO(6-5) ALMA observations 
(Garcia-Burillo et al. 2014) and which extends up to 2 arcsec (140 
pc) from the nucleus (see Fig. 5, bottom-right panel). 

Muller Sanchez et al. (2009) report the presence of a linear 
structure or northern tongue of H 2 emission streaming towards the 
nucleus. A good correlation between this northern tongue of H 2 
gas and dust emission at 12/im was presented by these authors 
(also seen in the dust maps of Prieto et al. 2014), which provides 
further evidence that the tongue is transporting feeding material 
to the AGN (Tomono et al. 2006; though see Barbosa et al. 2014 
for an outflow scenario). The superposition of the warm H 2 gas 
contours on the V — // colour dust map (Fig. 5, middle-left panel) 
shows that the warm molecular gas is not only coincident with dust 
emission in the northern tongue but that an overall clear correlation 
is observed between the dust morphology and the H 2 emission, 
including the gas concentrations east and west of the nucleus and 
the ring-like structure ~2 arcsec south of it. The correlation is 
however not one to one, as the warm H 2 gas presents a wider and 
more extended distribution than the dust. Despite the wider FoV of 
NaCo, no molecular H 2 gas emission is observed at scales larger 
than 2 arcsec, which indicates that all the warm molecular H 2 gas 
of NGC 1068 is concentrated in the central 140 pc. 

In Section 4.4 we found that the X-ray emission is the most 
likely cause of excitation of the nuclear warm H 2 molecular gas 
for most sources. In some regions close to the strong nuclear X-ray 
source and devoid of dust the X-ray emission could be able not 
only to excite the H 2 molecules but even to dissociate them. This 
extreme case might be observed for NGC 1068, showcasing the 
role of dust in shielding H 2 molecules from the intense X-ray 
radiation: the south-west cavity seen in the dust map and H 2 
contours (Fig. 5, middle-left panel) coincides with an excess of 
V-band (HST/F550W; Fig. 5, top-left panel) and Chandra emission 
(Fig. 5, middle-right panel), while the Chandra knot ~1 arcsec 
north of the nucleus coincides with a depression in the dust map 
and a hole in the H 2 molecular gas distribution (Fig. 5, top-right 
panel); finally, an arch of X-ray Chandra emission ~0.5 arcsec east 
of the nucleus is observed to coincide with an equivalent depression 
zone in the dust and H 2 gas maps. For NGC 1068 we might thus be 
seeing evidence that in the regions devoid of dust the H 2 molecule 
is not protected from the X-ray radiation but destroyed by it. This 
would explain the lack of H 2 emission in those regions where the 
X-ray emission is stronger, further supporting the existence of a 
clear link between H 2 molecular gas and X-ray radiation. 

NGC 1097 

NGC 1097 is a LINER/Syl galaxy located at 14.2 Mpc (Tully 
et al. 2009) with a prototypical nuclear ring (diameter ~18 
arcsec; ~ 1.2 kpc). AO IR broad-band imaging by Prieto et al. 
(2005) revealed dusty spirals around the nucleus and reaching the 
star-forming ring. 

The H 2 image of NGC 1097 (Fig. 6, top-right panel) is 
dominated by the nuclear emission, which can be traced up to 2 
arcsec (~ 140 pc) from the nucleus. Further out, in the star-forming 
ring at r ~ 8 arcsec, several compact emission regions are detected. 
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These H 2 point sources are coincident with the stellar clusters in the 
K-band (Fig. 6, top-right panel; see also Prieto et al. 2005). 

The nuclear H 2 emission shows a slightly elongated 
morphology along the east-west axis (middle-left panel in Fig. 6). 
This is in contrast with the H 2 1-0 S(l) intensity map obtained by 
Davies et al. (2009) with SINFONI (Fig. 6, middle-right panel), 
which shows a more compact morphology elongated along the 
north-south direction. The discrepancy is probably ascribed to the 
method used to estimate the continuum emission in Davies et al. 
(2009): the stellar and non-stellar continua are separated prior to 
the H 2 extraction, based on the dilution of the 2.3 \im CO band 
head. In order to accommodate the differences between both maps, 
the non-stellar emission should be also more extended along the 
north-south axis, but this is not observed neither in the Ks-brnd 
continuum nor in the I — K colour map. The morphology of dust 
as traced by ihd I —K colour (Fig. 6, middle-left panel) does in any 
case not correlate with the nuclear extent of the warm H 2 emission. 


NGC 1386 

NGC 1386 is a spiral Sb(s)a galaxy located in the Fornax cluster, 
at a distance of 15.3 Mpc (Jensen et al. 2003), that hosts a Sy2 
nucleus. NGC 1386 has a star-forming ring with a diameter of 
30 arcsec (2.2 kpc) and tightly collimated ionized gas in the HST 
[O III] and Ka images extending 2 arcsec-3 arcsec towards north 
and south from the optical nucleus (Ferruit et al. 2000; Prieto et al. 
2014). It shows a one-sided jet oriented in the south-east direction 
(Nagar & Wilson 1999; Fig. 7, middle-left panel). The nucleus (or 
IR peak; Prieto et al. 2014) is obscured, located 0.23 arcsec north 
of the optical peak of emission. 

The H 2 emission traced by NaCo presents a strong nuclear 
component and two extended regions up to 4 arcsec (~ 320 pc) 
towards the north and south of the nucleus (Fig. 7, top-right 
panel). The latter were reported by Reunanen et al. (2002) using 
long-slit spectroscopy in the near-IR. Its morphology, unveiled in 
this work for the first time, resembles that of the emission-line 
plumes detected in [O III] and [N ll]-i-Ha by Ferruit et al. (2000). 
In the / — K colour map (Fig. 7, middle-left panel), a dust lane 
is visible extending ~ 3 arcsec towards north-east and south-west 
of the infrared nucleus. The peak of the H 2 emission is associated 
with the nucleus. The morphology of the extended H 2 emission 
resembles that of the dust lanes covering the IR nucleus, however it 
presents a much wider extent (middle-left panel in Fig. 7). No H 2 
emission is detected in the star-forming ring. 


NGC 1566 

NGC 1566, at a distance of 20.5 Mpc (taken from NED adopting 
Ho = 73 km s“^ Mpc“^), is a nearly face-on SABbc spiral galaxy. 
It hosts a Seyfert nucleus that varies from type 1 to 2 (Alloin et al. 
1985) and has a one-sided ionization cone detected in [Olll] that 
extends 0.3 arcsec towards south-east (Schmitt & Kinney 1996); 
see also Fig. 8, middle-right panel. The optical HST images reveal 
dusty spirals reaching towards the nucleus (Malkan et al. 1998; see 
also dust map in Fig. 8, middle-left panel). Recent high-resolution 
CO ALMA observations show that the cold molecular CO(3-2) gas 
emission extends up to 3.5 arcsec forming a 2-arm nuclear spiral 
structure that corresponds very well with the dust lanes observed in 
the HST images (Combes et al. 2014) and the dust map in Fig. 8 
(bottom-left panel). 

The NaCo warm H 2 molecular gas emission of NGC 1566 


reveals an almost circular diffuse morphology (Fig. 8, top-right 
panel), slightly elongated in the east-west direction at r < 1 
arcsec, similar to the H 2 emission detected by SINFONI (Fig. 8, 
middle-right panel). The extent of the warm H 2 gas region in the 
NaCo image is 2 arcsec (~220 pc), which is ~2 times the size of 
that detected in the instrument-limited FoV of SINFONI and nearly 
the same as the extension of the H 2 emission detected by near-IR 
long-slit spectroscopy (Reunanen et al. 2002). The warm H 2 gas 
emission does not follow the dust spiral structure as does the CO 
emission (Fig. 8, bottom panels). 

NGC 7582 

NGC 7582 is a SBab galaxy located at 19.9 Mpc in the Grus 
Quartet (Terry et al. 2002). Its Sy2 nucleus, a strong IR source 
only seen from the /-band on (Prieto et al. 2014), is surrounded 
by star-forming clusters resolved in the nuclear starburst (Prieto 
et al. 2002; also visible in the VISIR/PAH2J2 image (11.88 pm) in 
Fig. 9). The morphology of the warm H 2 gas emission, traced by 
SINFONI, extends up to 2 arcsec (192 pc) and is associated with 
the starburst emission. The most embedded clusters, detected in the 
mid-IR and radio wavelengths (Wold & Galliano 2006; Orienti & 
Prieto 2010), show a counterpart in the H 2 map (Fig. 9 bottom-left 
panel). A compact unresolved source of warm H 2 gas is detected 
at the position of the AGN and is most likely excited by X-ray 
emission (Fig. 9, bottom-right panel). 


